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Dedicated to Richard Stanley on the occasion of his sixtieth birthday 

Abstract. It is proved that a certain symmetric sequence (ho, hi, ■ ■ ■ , hd) 
of nonnegative integers arising in the enumeration of magic squares of given 
size n by row sums or, equivalently, in the generating function of the Ehrhart 
polynomial of the polytope of doubly stochastic n X n matrices, is equal to 
the h-vector of a simplicial polytope and hence that it satisfies the conditions 

of the g-theorem. The unimodality of (ho, hi, , h^), which follows, was 

conjectured by Stanley (1983). Several generalizations are given. 



A magic square is a square matrix with nonnegative integer entries having all 
line sums equal to each other, where a line is a row or a column. Let H n (r) be the 
number ofnxn magic squares with line sums equal to r. The problem to deter- 
mine H n (r) appeared early in the twentieth century ^UJ- Since then it has attracted 
considerable attention within areas such as combinatorics, combinatorial and com- 
putational commutative algebra, discrete and computational geometry, probability 
and statistics [H EI IH El EI CE1 113 IIHl EH EH| It was conjectured by Anand, 
Dumir and Gupta and proved by Ehrhart |Hj and Stanley ^7] (see also |23l 
Section 1.5] and |241 Section 4.6]) that for any fixed positive integer n, the quantity 
H n (r) is a polynomial in r of degree (n — l) 2 . More precisely, the following theorem 
holds. 

Theorem 1.1. (Stanley [T?|[T§]) For any positive integer n we have 



where d = n 2 — 3n + 2 and the hi are nonnegative integers satisfying ho — 1 and 
hi = hd-i for all i. 

It is the first conjecture stated in [231 (see Section 1.1 there) that the integers hi 
appearing in satisfy further the inequalities 



In this paper we prove this conjecture by showing that (h , h\, . . . , hd) is equal to 
the h-vectov of a d-dimensional simplicial polytope. Such vectors are known to be 
symmetric and unimodal and are characterized by McMullen's ^-theorem see 
[1EH and Section 2. 
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(1) 




(2) 



ho < hi < • • • < h 
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A few comments on Stanley's conjecture and the proof given in this paper are in 
order. It is known that the sequence (ho, hi, ... , hd) of Theorem ll.ll is a Gorenstein 
sequence, meaning it is the /i-vector of a standard, graded, Gorenstein commutative 
ring; see for instance |231 Section 1.13]. It is an important open problem to charac- 
terize Gorenstein sequences; see |231 Section II. 6]. In this direction it was originally 
conjectured by Stanley |18j that a sequence (ho, hi, . . . ,hd) is Gorenstein if and 
only if it satisfies the conditions of the (/-theorem but a counterexample was later 
given in [20]. Our result ( Corollary 13. 6fl and its generalization to the enumeration 
of magic labelings of regular bipartite graphs (Corollary 14. 5f) give an instance in 
which Stanley's original conjecture turns out to be true. Another such instance, in 
which the entries of (ho, hi, ... , hd) count linear extensions of a naturally labeled 
poset by the number of descents, was given recently by Reiner and Welker [15] . 
The polynomial H n (r) is the Ehrhart polynomial (see Section^ of the Birkhoff 
polytope of doubly stochastic n x n matrices. We will state our main result in the 
context of Ehrhart polynomials of integer polytopes (Theorem 13. 50 as well as that 
of enumerating solutions to systems of linear homogeneous Diophantine equations 
(Corollary 14. 1[) and will show that both situations of Theorem 11.11 and [TB] appear 
as special cases. 

This paper was largely motivated by the work [15 of Reiner and Welker. I 
am grateful to Volkmar Welker for encouraging discussions and to Jesus DeLoera, 
Victor Reiner, Francisco Santos and Richard Stanley for helpful suggestions. 



In this section we review some basic definitions and background on convex poly- 
topes and their face numbers, triangulations and Ehrhart polynomials. We refer 
the reader to the texts by Stanley |23l [2] , Sturmfels [2Zj and Ziegler [3Hj for more 
information on these topics. We denote by N the set of nonnegative integers. 

Face enumeration. Given a finite (abstract or geometric) simplicial complex A 
of dimension d — 1, let fi denote the number of i-dimensional faces of A, so that 
(fo, fi, . . . , fd-i) is the f-vector of A. The polynomial 



where /_i — 1 unless A is empty, is the h- polynomial of A, denoted h(A,t). The 
h- vector of A is the sequence (ho, hi, ... , hd) defined by J^J. 

A polytopal complex T |28l Section 8.1] is a finite, nonempty collection of convex 
polytopes such that (i) any face of a polytope in T is also in T and (ii) the inter- 
section of any two polytopes in T is either empty or a face of both. The elements 
of T are its faces and those of dimension are its vertices. The dimension of T 
is the maximum dimension of a face. The complex T is pure if all maximal faces 
of T have the same dimension. The collection T(P) of all faces of a polytope P 
and the collection F(dP) of its proper faces are pure polytopal complexes called 
the face complex and boundary complex of P, respectively. Thus P is simplicial if 
J-(dP) is a simplicial complex. The /i-vectors of boundary complexes of simplicial 
polytopes are characterized by McMullen's g-theorem ^I] [231 Section III. 1] [2HI 
Section 8.6] as follows. A sequence (go, gi, ■ ■ ■ , gi) of nonnegative integers is said to 
be an M -vector if 
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(i) go — 1 and 

(ii) 0<g l+ i <gf } bal<i<i-l, 
where 0^ = and 

for the unique representation 

+ + - + (* 

with ki > ki—i > ■ ■ ■ > kj > j > 1 , if n > 1 . A sequence (Jig, hi, ■ ■ ■ , hd) of 
nonnegative integers is the /i-vector of the boundary complex of a d-dimensional 
simplicial polytope if and only if 

(i) hi = hd-i for all i and 

(ii) (ho, hi -h a ,..., fr|d/2J ~ h ld/2\-i) is an M- vector. 

In particular (ho, hi, . . . , hd) is symmetric and unimodal and hence satisfies the 
inequalities J3J, known as the Generalized Lower Bound Theorem for simplicial 
polytopes. 

THangulations and Ehrhart polynomials. A triangulation of a polytopal com- 
plex J 7 is a geometric simplicial complex A with vertices those of T and underlying 
space equal to the union of the faces of T ', such that every maximal face of A is 
contained in a face of T . A triangulation of the face complex T(P) of a polytope 
P is simply called a triangulation of P. 

For any set a consisting of vertices of the polytopal complex T we denote by 
T\a the subcomplex of faces of T which do not contain any of the vertices in a 
and write T\v for J-\a if a consists of a single vertex v. Given a linear ordering 
t = (vi,V2, ■ ■ ■ ,v p ) of the set of vertices of T we define the reverse lexicographic 
triangulation or pulling triangulation A (J 7 ) = A T (J^) with respect to r |22| [§| |27l 
p. 67] as A(T) = {v} if T consists of a single vertex v and 

A(^) = A(^\« p ) U (J {conv({v p } U G) : G E A(^(F)) U {0}} 

F 

otherwise, where the union runs through the facets F not containing v p of the 
maximal faces of J- which contain v p and A(jF\w p ) and A(JF(F)) are defined with 
respect to the linear orderings of the vertices of F\v p and F, respectively, induced 
by t. Equivalently, for io < ii < ■ ■ ■ < it the set {vi , Vi 1 , . . . , Vi t } is the vertex set 
of a maximal simplex of A r (J 7 ) if there exists a maximal flag Fq C Fi C • • • C Ft 
of faces of T such that Vi is the last vertex of Fj with respect to r for all j and Vi 
is not a vertex of -Fj-i for j > 1. A different way to define A T (J r ) is the following. 
For any vertex v of let 

pulUJF) = (^\«) U (J {conv({w} U G) : G e ^(F) U {0}}, 

F 

where the union runs through the facets -F not containing v of the maximal faces 
of T which contain v. If !Fq = T and T% — pull^ _ (Ti-i) for 1 < i < p then 
is a triangulation of T which coincides with A T (J-). It follows from |12l Theorem 
2.5.23] (see also p. 80]) that if T is the boundary complex of a polytope P then 
pull 1 ,(J r ) is the boundary complex of another polytope, obtained from P by moving 
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its vertex v beyond the hyperplanes supporting exactly those facets of P which 
contain v. This observation implies the following lemma. 

Lemma 2.1. The reverse lexicographic triangulation of the boundary complex of a 
polytope with respect to any ordering of its vertices is abstractly isomorphic to the 
boundary complex of a simplicial polytope of the same dimension. 

A convex polytope P C R 9 is said to be a rational or an integer polytope if 
all its vertices have rational or integer coordinates, respectively. It is called a 0-1 
polytope if all its vertices are 0-1 vectors in R 9 . If P is rational then the function 
defined for nonnegative integers r by the formula 

Ehrhart(P, r) = # (rP n Z 9 ) 

is a quasi-polynomial in r, called the Ehrhart quasi-polynomial of P |241 Section 
4.6]. If P is an integer polytope then this quasi-polynomial is actually a polynomial 
in r. Let A C R 9 be the affine span of the integer polytope P. A triangulation A of 
P is called unimodular if the vertex set of any maximal simplex of A is a basis of the 
affine integer lattice ^4nZ 9 . We denote by A T the reverse lexicographic triangulation 
of an arbitrary polytope P with respect to the ordering t of its vertices. Following 
[2"2] we call such an ordering of the vertices of an integer polytope P compressed 
if A T is unimodular and call P itself compressed if so is any linear ordering of 
its vertices. The following lemma holds for any unimodular triangulation of P, 
although we will not need this fact here. 

Lemma 2.2. ( [221 Corollary 2.5]) If P is an m- dimensional integer polytope in R 9 
and t is a compressed ordering of its vertices then 

Y Ehrhart(P,r)f = /^T'l, . 



If P C R m is an m-dimensional polytope and V is any linear subspace of R m 
then the quotient polytope P/V C R m / V is the image of P under the canonical 
surjection R m — > R m / V. This is a convex polytope in R m / V linearly isomorphic 
to the image 7r(P) of P under any linear surjection 7r : R m — > R m - dlmy with kernel 
V. Recall that the simplicial join Ai * A2 of two abstract simplicial complexes Ai 
and A2 on disjoint vertex sets has faces the sets of the form o\ U 0-2, where a\ € Ai 
and <72 G A2 and that h(A\ * A2,i) = h(Ai,i) /i(A2,t). The following proposition 
is essentially Proposition 3.12 in |15j . 

Proposition 2.3. Let P be an m-dimensional polytope in R m having a triangu- 
lation abstractly isomorphic to a * A, where a is the vertex set of a simplex not 
contained in the boundary of P. Let V be the linear subspace o/R m parallel to the 
affine span of a. 

The boundary complex of the quotient polytope P/V C R m / V is abstractly iso- 
morphic to J-{P)\a and inherits a triangulation abstractly isomorphic to A. 

Two compressed polytopes. (a) A real n x n matrix is said to be doubly sto- 
chastic if all its entries are nonnegative and all its rows and columns sum to 1. The 
set P of all real doubly stochastic n x n matrices is a convex polytope in R nxn of 
dimension (71 — l) 2 , called the Birkhoff polytope |28l Example 0.12]. It follows from 
the classical Birkhoff-von Neumann theorem that the vertices of P are the n x n 
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permutation matrices, so that P is a 0-1 polytope. The Birkhoff polytope was 
shown to be compressed by Stanley |22l Example 2.4 (b)] (see also |27l Corollary 
14.9]). 

(b) Let Q be a poset (short for partially ordered set) on the ground set [m] := 
{1,2,..., to}. Recall that an (order) ideal of is a subset I C ft for which i <n j 
and j s / imply that i e I. Let Q° be the poset obtained from f2 by adjoining a 
minimum element 0. The order polytope |25| of f2, denoted 0(Q), is the intersection 
of the hyperplane xq = 1 in K m+1 with the cone defined by the inequalities a;, > Xj 
for i < j in £7° and Xi > for all i. Thus O(fi) is an m-dimensional convex polytope. 
The vertices of O(fi) are the characteristic vectors of the nonempty ideals of H° 
[2"51 Corollary 1.3] so, in particular, 0(f2) is a 0-1 polytope (see [23 Theorem 1.2] 
for a complete description of the facial structure of 0(H)). Order polytopes were 
shown to be compresed by Ohsugi and Hibi p3 Example 1.3 (b)]. 

3. Special simplices 

Throughout this section P denotes an m-dimensional convex polytope in R 9 with 
face complex ^(P). Let E be a simplex spanned by n vertices of P. We call E a 
special simplex in P if each facet of P contains exactly n — 1 of the vertices of E. 
Note that, in particular, E is not contained in the boundary of P. 

Example 3.1. Let P be the polytope of real doubly stochastic n x n matrices. 
If V\, V2, ■ ■ ■ , v n are the n x n permutation matrices corresponding to the elements 
of the cyclic subgroup of the symmetric group generated by the cycle (12 • • • n) 
(or any n permutation matrices with pairwise disjoint supports) then Vi, v%, ■ ■ ■ , v n 
are the vertices of a special simplex in P. Indeed, each facet of P is defined by an 
equation of the form Xij = in K nxn and misses exactly one of Vi, V2, ■ ■ ■ , v n . 

Example 3.2. Let f2 be a poset on the ground set [m] := {1,2, ...,m} which 
is graded of rank n — 2 (we refer to [2] Chapter 3] for basic background and 
terminology on partially ordered sets) and P = 0(0,) be the order polytope of 
in ffi. m+1 . Let f2° be the poset obtained from O by adjoining a minimum element 
and for 1 < i < n let Vi be the characteristic vector of the ideal of elements of fl° 
of rank at most i — 1, so that Vi is a vertex of P. Since a facet of P is defined cither 
by an equation of the form Xi — Xj with i < j in fl° and i, j in successive ranks or 
by one of the form Xi = for i 6 il° of rank n — 1, it follows that V\, 1)2, ■ ■ ■ , v n are 
the vertices of a special simplex in P. 

Lemma 3.3. Suppose that v\, 1)2, ■ ■ ■ ,v n are the vertices of a special simplex in P. 
If F is a face of P of codimension k for some 1 < k < n — 1 and F does not contain 
any of V%, 1)2, ■ ■ ■ , Vk then F must contain Vi for all k + 1 < i < n. 

Proof. Let E be the special simplex with vertices v\, V2, ■ ■ ■ , v n . Any codimension 
k face of a polytope can be written as the intersection of k facets, so we can write 
F = F\ n F "2 n ■ ■ • n Fk where the Fj are facets of P. For each 1 < i < k we 
have Vi ^ F and hence ^ Fj for some j = ji. Since E is special the integers 
ji , j2 , . . . , jk are all distinct and hence for each 1 < j < k we have ^ Fj for some 
1 < i < k, which in turn implies that V{ S Fj for all k + 1 < i < n. It follows that 
Vi E Fx n F 2 n • • ■ n F k = F for all k + 1 < i < n. □ 

Lemma 3.4. Suppose that r = (v p , u p -i, . . . , v±) is an ordering of the vertices of P 
such that a = {v±, V2, ■ ■ ■ , v n } is the vertex set of a special simplex in P. Let A be 
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the abstract simplicial complex on {y n +i, ■ ■ ■ , v p } defined by the reverse lexicographic 
triangulation of J-(P)\a with respect to (v p ,v p -i, . . . ,v n +i)- 

(i) The reverse lexicographic triangulation A T of P is abstractly isomorphic to 
the simplicial join a * A. 

(ii) A is abstractly isomorphic to the boundary complex of a simplicial polytope 
of dimension m — n + 1. 

Proof, (i) Let &i — {vi, ...,Vi} for < i < n, so that ctq = and a n = a, and 
let A, denote the abstract simplicial complex on the set {"Ui+i, . . . ,v p } defined by 
the reverse lexicographic triangulation of T(P)\o~i with respect to the ordering 
{v p ,Vp-i, . . . , u<+i). To prove that A = a n * A n , which is the assertion in the 
lemma, we will prove that J-(P)\o~i is pure (m — i)-dimensional and that Ao = 
Oi * Aj for all < i < n by induction on i. This is obvious for i = so let 1 < i < n. 
By induction, any maximal face F of jF(P)\(T,_x is a codimension i — 1 face of 
P. Since F does not contain any of the vertices v\, . . . , Mt-i, by Lemma 13.31 we 
have Vi 6 F. This implies that ^F(P)\at is pure (m — «)-dimensional and that 
A.;_i = * Aj. The last equality and the induction hypothesis Ao = <Tj_i * Aj_i 
imply that Ao = o~i * Aj, which completes the induction. 

(ii) Let V be the linear subspace of parallel to the afhne span of the vertices in a 
and P/V be the corresponding quotient polytope of P, so that P/V has dimension 
m — n + 1. Part (i) and Proposition 12.31 imply that A is abstractly isomorphic 
to a reverse lexicographic triangulation of the boundary complex of P/V . This is 
in turn isomorphic to the boundary complex of a simplicial polytope of dimension 
m — n + 1 by Lemma 12.11 □ 

The following theorem is the key to the results in this paper. 

Theorem 3.5. Suppose that P is an integer polytope and r = (v p , v p -\ 1 . . . , v\) is 
an ordering of its vertices such that: 

(i) r is compressed and 

(ii) {^i, i>2, ■ ■ ■ , v n } is the vertex set of a special simplex in P. 



where h(t) = ho + h%t + • • ■ + hjt d is the h-polynomial of the boundary complex 

of a simplicial polytope Q of dimension d — m — n + 1, so that h(t) satisfies the 

conditions in the g-theorem. 

In particular hi = hd—i for all i and 1 = ho < hi < • ■ ■ < h^/2] ■ 

Moreover, Q can be chosen so that its boundary complex is abstractly isomorphic 

to the reverse lexicographic triangulation o/^ r (P)\ {v\, . . . , v n } with respect to the 

ordering (v p ,v p -i, . . . ,w„+i). 

Proof. Let a = {vi,V2, . . . ,f n } and let A denote the reverse lexicographic trian- 
gulation of J-(P) \ a with respect to the ordering (v p , v p -i, . . . , Lemma |2.2I 
guarantees that the proposed equation holds with hit) — ft.(A T ,i). Part (i) of 
Lemma T3.4I implies that 



since face complexes of simplices have /i-polynomial equal to 1, and the result 



Then 




h(A T , t) = h{a *A,t) = h(a, t) h(A, t) = h(A, t), 



follows from part (ii) of the same lemma. 



□ 
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We now apply Theorem 13.51 to the Birkhoff polytope and to order polytopes of 
graded posets. Observe that our theorem does not apply to all integer polytopes 
since 0-1 polytopes with no regular unimodular triangulations are known to exist 

El 

Magic squares and the Birkhoff polytope. Let P be the polytope of real 
doubly stochastic n x n matrices. Observe that the polynomial Ehrhart(P, r) co- 
incides with the function H n (r) of Theorem 11.11 Since P is a compressed integer 
polytope of dimension (n — l) 2 , Theorem 13.51 and Example 13.11 imply immediately 
the following corollary. 

Corollary 3.6. For any positive integer n we have 

y Hn {r)f = - — 

r>0 v ; 

where h(t) = ho + hit + • • ■ + hdt d is the h-polynomial of the boundary complex of a 
simplicial polytope of dimension d = n 2 — 3n + 2, so that h(t) satisfies the conditions 
in the g-theorem. 

In particular hi = h^-i for all i and 1 = ho < hi < • • • < fyd/aj ■ 

In view of the last statement in Corollary 14.11 the polytope in the previous 
corollary can be constructed by pulling in an arbitrary order the vertices of the 
quotient of P with respect to the affine span of the vertices V\ , v 2 , . . . , v n , chosen 
explicitly as in Example 13. II 

Eulerian polynomials and equatorial spheres. Let be a graded poset on the 
ground set [to] := {1, 2, . . . , to} of rank n — 2. Let fii be the set of elements of f2 of 
rank %— \ for 1 < i < n — 1 and £(fi) be the set of linear extensions of Q, meaning 
the set of permutations w — (w±,W2, ■ ■ ■ ,w m ) of [to] for which uii <o Wj implies 
i < j . We assume that £1 is naturally labeled, meaning that the identity permutation 
(1,2,..., to) is a linear extension. The Q-Eulerian polynomial is defined as 

w(n,t)= y tdcs{w) 

where 

des(iu) =#{ie[m-l] :iu, > Wi+i} 

is the number of descents of w. Following ^S] we call a function g : fl — > R equatorial 
if min a6 n g{p) = and for each 2 < i < n — 1 there exist Oi_i G and G Oj 
such that di-i <n a i and 9( a i-i) — g( a i)- An ideal / or, more generally, a strictly 
increasing chain of ideals I\ C I2 C • • • C Ik in f2 is equatorial if the characteristic 
function xi of / or the sum xii + XI2 + ' " ' + Xih 1 respectively, is equatorial. The 
equatorial complex A eq (£l), introduced in ^S], is the abstract simplicial complex 
on the vertex set of equatorial ideals of £1 whose simplices arc the equatorial chains 
of ideals in f2. 

The following theorem is proved in Corollary 3.8 and Theorem 3.14 of |15j . 

Theorem 3.7. ( Reiner- Welker j!5j ) Let Q be a naturally labeled, graded poset on 
[to] having n — 1 ranks. The equatorial complex A eg (f2) is abstractly isomorphic to 
the boundary complex of a simplicial polytope of dimension d = to — n + 1 which 
has h-polynomial equal to the fl-Eulerian polynomial W(Cl,t). 
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Hence W(fl, t) satisfies the conditions in the g-theorem and, in particular, it has 
symmetric and unimodal coefficients. 

Let P be the order polytope of SI and Sl° be the poset obtained from SI by ad- 
joining a minimum clement 0. Recall that the vertices of P are the characteristic 
vectors of the nonempty ideals of Sl°. The order polytope comes with its canonical 
triangulation |25| [151 Proposition 2.1], which is a unimodular triangulation with 
maximal simplices bijecting to the linear extensions of SI. This canonical triangu- 
lation is in fact the reverse lexicographic triangulation of O(Sl) with respect to any 
ordering (u p) u p -\, . . . , uf) of its vertices such that i < j whenever the ideal of fi° 
defined by m is strictly contained in that defined by uj . We will use the following 
lemma. 

Lemma 3.8. Let Vj, be the characteristic vector of the ideal of elements of 0° of 
rank at most i—1 for 1 < i < n. Let a — {v\, . . . , v n } and r = (v p , . . . , be an 

ordering of the remaining vertices of P such that i < j whenever i, j • > n + 1 and 
the ideal defined by Vi is strictly contained in that defined by Vj . 

The equatorial complex A eg (Sl) is the abstract simplicial complex defined by the 
reverse lexicographic triangulation of J-(P)\o~ with respect to r. 

Proof. Let T denote the face complex of P and let xj = by convention. The 
maximal faces of F\a are the faces of P defined by systems of equations of the 
form Xi s — x.j 3 for < s < n — 1 where (i) «o = and j n -i = L (ii) *s £ for 
1 < s < n — 1, j s € r^s+i for < s < n — 2 and i s <q j s for 1 < s < n — 2 and 
(iii) if j s = i s +i for consecutive values s = a, a + 1, . . . , b — 1 of s then the interval 
[i a ,jb] in SI consists only of the elements of the chain i a < i a +\ < ■ • • < ib < jb- 
The statement of the Lemma follows from the description of the maximal faces of 
a reverse lexicographic triangulation A(JF) (see Section|2Jl and that of the maximal 
faces of A eg (f2) (see ^1 Proposition 3.5]). We omit the details. □ 

Proof of Theorem \3.1\ Let P be the order polytope of fi, as before. Observe that 
Ehrhart(P, r) is equal to the number of order reversing maps p : O — > {0, 1, . . . , r}. 
It follows from j2U Theorem 4.5.14] that 



Let the vertices Vi, V2, ■ ■ ■ , v p of P and r = (v p , . . . , be as in Lemma 13.81 We 

checked in Example 13.21 that v% , v%, . . . , v n are the vertices of a special simplex in 
P. Since P is a compressed integer polytope (see Section |2J) Theorem 13.51 applies 
and we have 



where hit) = ho + hit + • • ■ + hdt d is the /i-polynomial of a simplicial polytope 
of dimension d = m — n + 1 having, in view of Lemma 13.81 boundary complex 
abstractly isomorphic to A eg (f2). Comparison with yields h(f) — W(fl,t) and 





completes the proof. 



□ 
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4. Rational polyhedral cones 



In this section we state several corollaries of Theorem 13.51 including a general- 
ization of Corollary 13.61 to magic labelings of bipartite graphs. Let $ be a p x q 
integer matrix with rank p and P$ be the monoid of vectors x € N 9 satisfying 
the homogeneous system of linear equations Qx = 0. We assume that £$ is 
nonzero, so that q > p + 1, and set m = q — p — 1. Wc denote by C$ the cone 
of vectors x — (x\,X2, ■ ■ ■ ,x q ) € R 9 satisfying <j>x = and Xi > for all i. Let 
L(x) = aixi + C12X2 + ■ • • + a q x q be a linear functional on R 9 with a, £ Q for all i 
such that the set 



is nonempty and bounded. Thus P is an m-dimensional convex polytope. Observe 
that 



for any nonnegative integer r. Let E$ be the submonoid of E$ consisting of those 
elements having positive coordinates. We say that (3 6 P$ is the unique minimal 
element of P$ if we have (3 < 7 coordinatewise for all 7 6 £"$ . 

Corollary 4.1. Let P &e an m-dimensional polytope, as in l|5*|l. J/ 

(i) P is an integer polytope, 

(ii) P$ /las a unique minimal element (3 and 

(hi) £/iere exists a compressed ordering t = (v p , v p -\, . . . , v\) of the vertices of 
P such that v\ + V2 + ■ ■ ■ + v n — [3 for some n 

then the conclusion of Theorem Iff. 51 holds. 

Proof. In view of Theorem 13. 51 it suffices to show that {v±, V2, ■ • ■ , v n } is the vertex 
set of a special simplex in P. Let (3 — (j3i, 02, ■ ■ ■ , P q ) and let F be a facet of P, 
so that F is defined by an equation of the form Xk = for some 1 < k < q. Wc 
need to show that exactly one of v±, V2, ■ ■ ■ ,v n has positive feth coordinate. Clearly 
at least one of v\, 1*2 , • ■ ■ , v n has this property, since Pi > for all i. Assume on 
the contrary that at least two Vj have positive fcth coordinate, say v\ and V2, so 
that 1 < 7fe < (3k if «i = (71, 72, ■ ■ ■ , 7 9 )- Since F is a facet of P there exists 
a point x = (x\, X2, ■ ■ ■ , x q ) in the affine span of P, which we may assume to be 
rational, satisfying Xk < and Xi > for all % ^ k. By replacing x with a suitable 
integer multiple we find a point a — (ai, CK2, • • • , ctq) € satisfying $ce = 0, 
afe < and a, > for all i ^ k. We may choose a nonnegative integer £ so that 
< a k + (3k + £7fc < (3k (with t = if a k + (3k > 0). Then a + /3 + tv\ is in P$ and 
has fcth coordinate strictly less than (3k, which contradicts the minimality of (3. □ 

The condition that P$ has a unique minimal clement is satisfied if (1, 1, . . . , 1) € 
R 9 is in P$ and is known to hold if and only if the semigroup ring P$ = k[P$], 
generated over a field k by the monomials corresponding to elements of P$, is 
Gorenstein; see for instance |23l Section 1.13]. 

Corollary 4.2. Let P be an m-dimensional polytope, as in JSJ. If 

(i) P is a compressed integer polytope, 

(ii) P$ is Gorenstein and 

(iii) P$ is generated as a monoid by the vertices of P 



(5) 



P = {x E C$ : L(x) = 1} 



(G) 



Ehrhart(P, r) = # {x 6 P$ : L(a;) = r} 
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then the conclusion of Theorem holds where, in the statement of the theorem, 
n has the value L(f3) for the unique minimal element [3 of E$. 

Proof. Let (3 be the unique minimal element of £$, whose existence is guaranteed 
by (ii) . Then (iii) implies that (3 = v\ + «2 + • ■ ■ + v n for some vertices v\ , i>2, ■ ■ ■ , v n 
of P, which must be pairwise distinct. Because of (i) any ordering (v p , v p -i, . . . , Vx) 
of the vertices of P satisfies the assumptions of Corollary 14. II The result follows 
from this corollary observing that L{(3) = n. □ 

General conditions on <I> and P which guarantee assumptions (i) and (iii) of 
Corollary 14.21 were given by Ohsugi and Hibi j!3| . 

Corollary 4.3. Let P be an m- dimensional polytope, as in If P is a 0-1 

polytope and i?$ is Gorenstein then the conclusion of Corollary \4-^\ holds. 



Proof. Since P is a 0-1 polytope, Theorem 1.1 and Lemma 2.1 in |13j imply, respec- 
tively, that P is compressed and that E$ is generated as a monoid by the integer 
points of P. which are exactly the vertices of P. Thus the assumptions of Corollary 
IP hold. □ 

Remark 4.4. Recall that a matrix is called totally unimodular if all its subdc- 
tcrminants are equal to 0,-1 or 1. Suppose that the linear functional L(x) — 
a\x\ + a2X2 + • • • + a q Xq satisfies a« 6 Z for all i. It follows from the results of 
[7| that P is an integer polytope if the matrix obtained from $ by adding the row 
(di, 0,2, . . . , a q ) is totally unimodular. Hence, in view of Corollary 14.31 this state- 
ment and the assumptions that (i) P C [0, l] q and (ii) i?$ is Gorenstein imply the 
conclusion of Corollary 14. II 



Magic labelings of graphs. Let G be a graph (multiple edges and loops allowed) 
with p vertices and q edges and edge set £ . A magic labeling ^7] of G of index r 
is an assignment £ : £ — > N of nonnegative integers to the edges of G such that for 
each vertex v of G the sum of the labels of all edges incident to v is equal to r, in 
other words, 

e: v£e 

Let M denote the real vector space with basis £ and let x e be the linear functional 
on M £ dual to the basis element e E £. Let $ be a full rank, p' x q integer matrix 
with kernel the set of x € Mr satisfying the linear system of equations of the form 

(7) E x * = E Xe ' 

e: v£e e': v'£e f 

in M. £ , where v, v' are vertices of G. If L{x) is any of the functionals on IR £ in @ and 
P is as in JSJ) then Ehrhart(P, r) counts the number Hc(r) of magic labelings of G of 
index r. It follows from |17l Proposition 2.9] and either 22, Theorem 2.3] (applied 
as in |22l Example 2.4 (b)] in the case of the Birkhoff polytope) or |13l Theorem 
1.1] that conditions (i) and (iii) of Corollary 14.21 are both satisfied if the graph 
G is bipartite (or, more generally, if it satisfies condition (iii) of |17l Proposition 
2.9]). Clearly condition (ii) is satisfied if G is regular, since then (1, 1, .... 1) e E<j>. 
Assuming further that G is connected we have p' — p — 2. The following corollary 
specializes to Corollary 13 . 61 when G is the complete bipartite graph on two sets of 
vertices, each of size n. 



MAGIC SQUARES AND A CONJECTURE OF STANLEY 
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Corollary 4.5. For n > 1 and for any connected regular bipartite graph G with p 
vertices and q = np/2 edges we have 

YH G (r)t r = - m 

r>0 V ; 

where m = q — p + 1 and = ho + hit + • • • + hdt d is the h-polynomial of the 
boundary complex of a simplicial polytope of dimension d — m — n + 1, so that h(t) 
satisfies the conditions in the g-theorem. 

In particular hi = h^-i for all i and 1 = ho < h\ < ■ ■ ■ < hy d / 2 \ ■ 
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